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The curl modes of cosmic microwave background polarization allow one to indirectly constrain 
the primordial background of gravitational waves with wavelengths roughly the horizon size or 
larger with frequencies below 10 -16 Hz. The planned high precision timing observations of a large 
sample of millisecond pulsars with the Pulsar Timing Array or with the Square Kilometer Array 
can either detect or constrain the stochastic gravitational wave background at frequencies greater 
than roughly 0.1 years^ 1 . While there are no strong observational constraints on the gravitational 
wave background across six or more orders of magnitude between 10 -16 Hz and 10 -10 Hz and 
it is difficult to get a constraint below 10 -12 Hz using objects in our Galaxy, we suggest that the 
anisotropy pattern of time variation of the redshift related to a sample of high redshift objects can be 
used to constrain the gravitational wave background around 10 -12 Hz. Useful observations for the 
monitoring of an anisotropy signal in a global redshift change include spectroscopic observations of 
the Ly-a forest in absorption towards a sample of quasars, redshifted 21 cm line observations either 
in absorption or emission towards a sample of neutral HI regions before or during reionization, and 
high frequency (0.1 Hz to 1 Hz) gravitational wave analysis of a sample of neutron star — neutron 
star binaries detected with gravitational wave instruments such as the Decihertz Interferometer 
Gravitational Wave Observatory (DECIGO). The low frequency background can also be constrained 
by arcsecond-scale anisotropy observations of the CMB. For reasonable observations in the future 
involving extragalactic sources, we find best limits at the level of Ugw < 10 _J at a frequency around 
10~ 12 Hz while the eventual ultimate limit one cannot beat is Ogw < 10 -11 . 

PACS numbers: 95.85.Sz 04.80.Nn, 97.10.Vm 



I. INTRODUCTION 

The observation of cosmic microwave background 
(CMB) anisotropies, especially using a map of the curl 
modes of polarization, allows one to constraint the pri- 
mordial background of gravitational waves with wave- 
lengths roughly the horizon size or larger with frequen- 
cies below 10 -16 Hz Recent temperature anisotropy 
observations at large angular scales with Wilkinson Mi- 
crowave Anisotropy Probe (WMAP; 0) yields a limit 
(at the 2cr level) of tt GW < 5 x 1CT 11 HQ with a maxi- 
mal sensitivity around 10~ 17 Hz. Here, f^cw is the frac- 
tional density contribution from a background of stochas- 
tic gravitational waves with the density pgw : 

M/) = -^, (i) 

Pc a log/ 

when the closure density is p c — 3Hq/8ttG. Here, Hq 
is the Hubble constant. At these super-horizon and 
horizon-size scales, a stochastic background of gravita- 
tional waves is expected from inflationary physics 0, |(| . 
The direct detection of this signal is now considered to 
be one of the primary goals of upcoming CMB polariza- 
tion anisotropy observations both from ground and space 
with planned missions such as the Inflation Probe of the 
NASA's Beyond Einstein Foundation Science Program. 
Relative to the current limit from intensity anisotropies, 
one can improve by at least 5 orders of magnitude to 
the ultimate level allowed by the cosmic variance-limited 
removal of a confusion related to lensing conversion of 
polarization modes Q • With a separation of lensing pro- 
duced curl-modes, we expect the limit to improve down 



to ficw < 10~ 16 with the best sensitivity again around 
a frequency around 10~ 17 Hz |8|. The proposed Explorer 
Probe of Inflationary Cosmology (EPIC), to be launched 
around 2012 as part of NASA's Beyond Einstein Foun- 
dation Science program, will be optimized to reach this 
level. 

When considering a subhorizon-scale wavelengths, the 
planned high precision timing observations of a large 
sample of millisecond pulsars with the Pulsar Timing 
Array or with the Square Kilometer Array can either de- 
tect or constrain the stochastic gravitational wave back- 
ground at frequencies greater than roughly 1 /T" bs where 
T Q bs is the total observational duration of the pulsar sam- 
ple [9j. While the waves with wavelengths lower than 
1/Tobs linearly change the observed interval of the pulse, 
the lower frequency limit of 1/T b s for the pulsar timing 
method comes from the fact that one cannot distinguish 
their effect from the long-term intrinsic spin evolution 
of pulsars. Since observations are likely to be restricted 
to less than 100 years, the pulsar timing arrays cannot 
constrain the gravitational wave background below 10~ 10 
Hz. Observations of two millisecond pulsars over 17 years 
limit the background to be f^w^ 2 < 2 x 10~ 9 (at the 
95% confidence level) [13, EH, while with the Pulsar Tim- 
ing Array, using a larger sample of pulsars, the limit is 
expected to reach the level of 2 x 10 -13 9]. At these 
frequencies, the stochastic background of gravitational 
waves is expected to be primarily dominated by massive 
black hole binaries ^3 and either a detection or a tight 
constrain on the background is useful to understand the 
extent to which massive black holes merge at redshifts 
around unity. 
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Below frequencies probed by the millisecond pulsar 
timing method, a background of stochastic waves is ex- 
pected from models related to global cosmic strings, 
phase transitions in the early universe, and cosmic turbu- 
lence . Note also that the primordial background be- 
low ~ 10~ 10 Hz cannot be constrained by the abundance 
of the light elements predicted by big bang nucleosyn- 
thesis 14]. At present, there are no strong observational 
constraints on the gravitational wave background across 
six or more orders of magnitude between 10~ 16 Hz and 
1CP 10 Hz. The binary pulsars give a limit f^w^ 2 < 0.04 
for 10- n Hz < / < 4.4 x lCT 9 Hz, and n GW h 2 < 0.5 
for 10~ 12 Hz < f < 10~ n Hz from the time variation of 
their orbital period compared with prediction of general 
relativity |ll| . In future they might give a limit at the 
^Gwh 2 ~ 10 -4 level that is determined by uncertain- 
ties in the local acceleration of these binaries ^^]. Us- 
ing Galactic objects, it is difficult to get information on 
the background below / < 10~ 12 Hz as this corresponds 
roughly to the typical distance to galactic objects. The 
longer waves act on the target object and the observer 
in the same manner, and the measurement sensitivity 
for these low frequency waves is strongly suppressed due 
to a cancellation effect 15]. Therefore, it is essential 
that we use extra-Galactic objects to set a limit for the 
low frequency background below 10~ 12 Hz. The current 
limit comes from observations of quasar proper motions 
across the sky (transverse motions) |l6|. For data col- 
lected over a 10 year span, the resulting constraint on 
the stochastic background at a frequency below 2 x 10~ 9 
Hz is ficw^ 2 < 0.11 ^3 ( see a ^ so EH f° r prospects with 
the Square Kilometer Array). 

In this paper, to probe the background at frequencies 
below 10 -12 Hz, we propose measuring the anisotropy 
pattern of time variation of the redshift related to a sam- 
ple of high redshift objects. Previously, observations of 
the time variation of the redshift have been proposed to 
study global cosmological parameters [l^. While cos- 
mology can induce either an acceleration or a deceler- 
ation to the redshift of an object, the resulting change 
can be extracted through the monopole change from a 
sample of objects spread over the sky since all sources 
are equally affected by cosmology. The resulting change 
due to a gravitational wave background, at the lowest 
order in spherical moments, is present in the quadrupole 
moment of the time evolving redshift [2(j. The use of 
the anisotropy pattern also allows one to separate sec- 
ondary effects that can cause a time-dependent redshift, 
such as the peculiar acceleration of the observer which is 
present via a dipole at the lowest order. We focus on the 
quadrupole measurement, but the measurement of addi- 
tional multiple moments of the anisotropy can be used 
to increase the confidence level of the derived limit. As 
we discussed, the proposed method with using the time 
variation of redshifts would be performed as a byprod- 
uct of the dark energy study, in contrast to the method 
based on the transverse motions that is basically caused 
by perturbative effects. Therefore we will pay special at- 



tention to relate the sensitivities of the measurement for 
the dark energy and the gravitational wave background. 
As we are studying cosmology, we can expect, in order 
of magnitude sense, that there is a simple relation be- 
tween their sensitivities that are properly normalized by 
cosmological parameters such as Hq. It is, however, im- 
portant to calculate the relevant coefficients in discussing 
future prospects. In fact, we find that these coefficients 
are largely different from unity. 

Here, extending the discussion in Ref. we con- 

sider a measurement of the quadrupole related to the 
anisotropy of the time derivative of Ly-a forest redshifts 
seen in absorption towards a sample of quasars. Addi- 
tional probes for the same purpose include low-frequency 
redshifted-21 cm observations either in absorption or 
emission towards a sample of neutral HI regions before 
or during reionization 21]. Another possibility is the 
high frequency (0.1 Hz to 1 Hz) gravitational wave anal- 
ysis of of a sample of neutron star-neutron star bina- 
ries detected with gravitational wave instruments such 
as the Decihertz Interferometer Gravitational Wave Ob- 
servatory (DECIGO) |22( . These observations have been 
suggested for measurement of cosmological parameters 
based on the phase variation induced during the propa- 
gation of gravitational waves in a cosmologically expand- 
ing background [23, |23| . While cosmologically induces 
a global phase shift, a low frequency gravitational wave 
background induces an anisotropy to this phase shift and 
its presence can be extracted from a sample of binaries 
spread over the sky with well- modeled wave forms. For 
reasonable future observations along these lines, we find 
best constraints at the level of Hqw < 10~ 5 at a fre- 
quency around 10~ 12 Hz, though, the eventual limit is 
probably at Slew < 10~ . Reaching this level is chal- 
lenging as it requires high resolution observations of 21 
cm background to resolve all halos spatially with fol- 
lowup high resolution spectral measurements over a 10 
year span over a bandwidth less than a kHz. In figure 
1, we plot sensitivities of various current and future con- 
straints for the low frequency gravitational wave back- 
ground. From this figure, the importance of the extra- 
Galactic probes to study waves below 10~ 12 Hz, but 
above the regime studied by CMB experiments, is clear. 

Incidently, one can also, in principle, use the monopole 
to constrain the energy density of gravitational back- 
ground. The resulting correction comes from the fact 
that a subhorizon-scale gravitational wave background 
acts as an additional radiation field in the universe and 
the expansion rate of the universe is modified to become 

h 2 (z) = Hq [o ro (i + zf + (n R + n GW )(i + z f + n A ] , 

for a spatially-flat universe with cosmological constant 
JIa, the matter density parameter fij^f and the radia- 
tion content flu, all relative to the critical density. Since 
matter and dark energy densities dominate the expansion 
rate over the redshift ranges one can probe with current 
observational techniques, we do not pursue a limit based 
on this argument further. An interesting constraint, how- 
ever, could potentially be derived if one were to study 
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FIG. 1: Sensitivities for various current (thin lines) 
and future (thick lines) measurements. Curves labeled 
"COBE/WMAP" and "EPIC/CMBpol" show the current 
and future constraints based on CMB anisotropics at large 
angular scales, while the limit labeled "CMB-small" is a po- 
tential constraint based on arcsecond resolution, or multipole 
moments between 5 x 10 4 to 5 x 10 6 , observations of CMB 
anisotropies with a residual radio point source contamination 
consistent with current anisotropy data. The curve labeled 
"binary psr" is the limit based on a sample of binary pulsars, 
based on the time variation of their orbits, the curve labeled 
"msp" is the current limit based on 17-year timing monitoring 
of 2 millisecond pulsars while the curve labeled "PTA" is the 
limit based on proposed Pulsar Timing Array network. The 
limit based on "QSO astrometry" involves proper motion ob- 
servations of a sample of radio bright quasars while "redshift 
variation" involves time monitoring of redshift variations in 
the Lyman-a forest towards bright quasars, the 21 cm line 
at very low radio frequencies and cosmological neutron star 
binaries with ultimate DECIGO. The curve labeled "21cm- 
ultimate" is the eventual ultimate limit for method with 21 
cm line. 



earliest times in the universe such as the era of big bang 
nucleosynthesis and using the number of massless neu- 
trinos. The constraint applies to waves with wavelengths 
below about a 10 pc, and the previous limit based on this 
argument constraints the high-frequency background, as 
observed today, to be f!cw < 5 x 10~ 4 with a matter- 
radiation equality at a redshift of ~ 3500 3J. With im- 
proving measurements on the ratio of photon-to-baryons 
from both the cosmic microwave background anisotropy 
spectrum, with physics dominating at the last scatter- 
ing surface around z ~ 1100, and the same ratio mea- 
sured independently based on nucleosynthesis arguments, 



and can derive an improved constraint on the subhori- 
zon gravitational background that existed between an age 
for the Universe of a minute to a few hundred thousand 
years. We will return to this topic in an upcoming paper. 

Just as horizon-size waves produce CMB anisotropies, 
a gravitational wave background at frequencies between 
10~ 14 to 10~ 12 Hz can produce fluctuations in the CMB 
temperature at arcsecond angular scales. Currently 
most CMB anisotropy observations probe down to ar- 
cminute scales from acoustic peak structure at the de- 
gree scale, and prospects for even further small-scale 
CMB anisotropies observations are limited. These ob- 
servations are also likely to be challenging given addi- 
tional secondary effects such as the Sunyaev-Zel'dovich 
effect |24j. While most secondaries can be separated, 
through properties such as the frequency spectrum rel- 
ative to CMB |25j, the foreground contamination from 
unresolved point-sources is expected to be the biggest 
confusion. We make an estimate on the expected con- 
straint by extrapolating the foreground contamination 
seen in current anisotropy data at arcminute scales. 

The discussion is organized as follows: In the next Sec- 
tion, we outline the calculation related to an anisotropy 
pattern in the time evolving redshift of a sample of 
sources in the presence of a gravitational wave back- 
ground. In § III, we discuss various cosmological probes 
related to this effect and estimate resulting limits that 
one can potentially obtain in the future related to the 
presence of a low frequency, but sub-horizon, background 
of gravitational waves. We conclude with a summary in 
§ IV. 

II. FORMULATION 

We first study how a background of stochastic gravita- 
tional waves affect the redshift of a photon, or a massless 
particle, propagating to an observer today from a source 
at a high redshift. We consider a spatially- flat Universe 
and write the metric that is perturbed by gravitational 
waves as 

ds 2 = a(r]) 2 {-drj 2 + [<% + h tj {-q, x)]dx i dx j } , (2) 

where r\ is the conformal time and the scale factor a(rj) is 
normalized such that 0(770) = 1 at the present epoch with 
r\ = 7/0. We write the conformal time interval and the 
normal time interval today as A77 and AT, respectively. 
Comparing to a background universe with hij = 0, the 
propagation time of a photon changes by 5rj due to the 
presence of gravitational waves. We can formally write 
down Sr/ as a line integral along the photon path: 

run- f Vo 

s V = -7r hijlv, {V ~ no)n]dr), (3) 

where n is the unit directional vector to the high redshift 
source and L is its comoving distance from the observer. 
This shift Sr] will result in an apparent time variation, 
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Az, to the redshift of the source, z, in a time interval 
AT: 



The angular power spectrum of coefficient zi m Gw is 

6Qgw 



Az 
AT 



» 9^ 
-i 9?7 



[»?, (V ~ r] Q )n}drj. 



(4) 



(2Zm;GW22'm';Gw) 



ir(l + 2)(l + 1)1(1-1) 



Since we are interested in wave lengths (< lOMpc) 
smaller than the horizon scale, we can safely express the 
gravitational waves for our study as 

hij(f,P,r),x) =exp[i27rf(r]-p-x)]D(T],f)bi j (fp), (5) 

where p is the unit vector for the propagation direction of 
the wave, D(i], /) represents the cosmological evolution 
of the wave amplitude and bij(fp) is the random number 
for each wave characterized by the vector fp. 
We decompose eq.JHJ as 

hij(f,p,r),x) = exp[j2vr/(?7 - p ■ x)} {D(r),f) - D(r] Q ,f)} 
+ exp[i2wf(r]-p-x)]D(r)o,f)b ij (fp) , (6) 

and substitute this in Eq. (0J • Our basic aim is to extract 
effects that are common to many sources. Here, the con- 
tribution of the first term in Eq. © is negligible since (i) 
we are considering a wave with a wavelength shorter than 
the horizon scale, (ii) the common coherent effect is made 
by the perturbation close to us, or the observer, and (iii) 
the evolution time scale of the amplitude (D/D)^ 1 is or- 
der of the age of the Universe. As is well known [2t| . 
the second term in Eq. JSJ can be expressed as a sum of 
two quantities, the information of the gravitational wave 
at the source and at the observer. The former can be 
dropped for the same reason as in the previous case. In 
summary, we can apply the standard argument for the 
low frequency gravitational background by using objects 
within our galaxy (e.g. timing residual for a galactic mil- 
lisecond pulsar) . Random variations in Az / AT for each 
binary due to the gravitational wave background can be 
regarded as a noise for measurement of the common sig- 
nal. For these binaries, however, the expected magnitude 
of variation is smaller than other sources of noise (e.g. lo- 
cal acceleration, detector noise) , and hereafter we neglect 
this noise component. 

Following the calculation of Burke [2(j , we obtain the 
angular pattern of the common redshift variation Az due 
to the background in the time interval AT as 



.(9) 

and the parameter flew here is defined by the normalized 
energy density Slow (/) per unit log frequency interval as 
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faut 



SlGw(f) 



(10) 



We have explicitly included a lower frequency cut-off 
since, below a certain frequency, energy density of the 
gravitational wave background will be severely con- 
strained by CMB polarization observations. 

We also define the angular power spectrum Ci that is 
useful for a statistical study as 



C, = 
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n(l + 2)(l + 1)1(1-1) 



(11) 



In eq.Q it is important to note that the effect of the 
gravitational wave background starts form quadrupole 
I = 2 mode and independent on the source redshift z, 
while the redshift change Az due to the cosmic accel- 
eration/deceleration is a monopole (I — 0) effect and 
depends on the redshift z 0] that corresponds to the 
coefficient ioo as 

Z00;acc = 

{ [q m (i + z) + n R (i + z f + n A (i + z)- 2 } 1/2 - 1} 

We define a parameter X that characterizes the second 
order correction to the relation between time intervals at 
the observer, At, and the source, At z due to effective 
acceleration as 



(1 + z)At z = At- XAt 2 + 0(At 3 



(12) 



The parameter X can have a finite value due to the cos- 
mological acceleration, gravitational wave background or 
the local acceleration. For each object we fit the param- 
eter X to statistically extract various information. For 
the cosmological acceleration we have 
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(13) 



where zgw(&t 4>), given through Eq. has the angular 
dependence of (1— cos 9) cos 2(f). This angular dependence 
can be written as a combination of Y\% and Y/_2, when 
I > 2 such that 



(1 — cos#)cos2</> 




(8) 



when the Hubble parameter at redshift z is H(z). The 
parameter X(z) can also be written as X(z) = 0.5[d(0) — 
a(z)], where the overdot represents a derivative with 
respect to the proper time. At small redshifts, when 
z « 1, X(z) — —Q.bqoz where go is the global deceler- 
ation parameter, Q m — fl\/2. In addition to time vari- 
ation in the redshift, for gravitational waves emitted by 
the source, the resulting effect related to global dynam- 
ics is an additional change to the phase of the Fourier 
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transformed gravitational waveform of a chirping binary 
by an amount oc f~ 13 / 3 X whose frequency dependence is 
largely different from the intrinsic binary evolution pre- 
dicted by the post Newtonian expansion. This phase 
correction has been used as a probe of cosmological pa- 
rameters in Ref. jl^, • The anisotropy of the phase 
correction, or anisotropy measurement of X from a large 
sample of binaries spread over the sky, can be used a 
probe of large wavelength gravitational waves. 

We assume that the number of the probes (e.g neutron 
star binaries, QSO line-of-sights) is sufficient and the ef- 
fects caused by a finite number of such sources, such as 
Poisson fluctuations, is not important for measurement of 
zim at low is. We also introduce the ratio R = X(z)/Hq 
to relate the resolution of the acceleration measurement, 
from which X(z) is derived, and the resulting constraint 
on the low-frequency gravitational wave background. In 
the case of the measurement related to the cosmic acceler- 
ation, in terms of the monopole variation associated with 
X(z), the redshift dependence of R is important and we 
analyze probes binned over some finite width in redshift. 
From the same data, but using higher order anisotropies, 
we get the constraint for FIgw- With a statistical re- 
lation between the measurement error of the spherical 
harmonic coefficients, zi m , in terms of an overall error in 
the power spectrum, we have 




where a 1 zoo is the variance related to the monopole mea- 
surement. The minimum energy density of the gravita- 
tional wave background measurable by any of the meth- 
ods discussed below can be obtained with 




To obtain the limit, under the null hypothesis of no gravi- 
tational wave background, we set C; = in Eq. [21 Since 
the measurement of the whole power spectrum related 
to redshift derivatives, or X(z), may not be achievable 
in the near future, we only concentrate on the variance 
related to the quadrupole, Ci . Note that the uncertainty 
in this measurement is now determined by the error to 
which the monopole can be established, with a 2 ZQ Qj given 
by the error in R = X(z)/H . In figure 2, we plot the 
angular spectrum C; for Flaw ■ We also put the measure- 
ment error for each C/ , if the sensitivity for the monopole 
mode R is 10~ 3 for a given redshift shell. 

Making use of the quadrupole anisotropy only, we ob- 
tain the 1-cr constraint for the background as 

Q GW < 4.0 x U) 2 (AR) 2 shell (16) 

from probes at each redshift shell, where AR is the er- 
ror to which R(z) can be established for such a shell. 
When we use the information of anisotropies up to I = 5, 
the coefficient 4.0 x 10 2 is reduced to 3.9 x 10 2 . As we 
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FIG. 2: The angular power spectrum Ci of a redshift shell 
due to the low frequency gravitational wave background with 
SJgw = 10 . For the error bars we assume a case that we 
can measure the monopole mode R with sensitivity 10 -3 for 
the shell. 

commented earlier, the signal due to the gravitational 
wave background does not depend on the redshift. This 
means that the actual constraint is better than the result 
using a single shell as one can combine shells at different 
redshifts, independent of the exact redshift. For cosmo- 
logical model with h — 0.7, fi m = 0.3, and Ha — 0.7 the 
magnitude of R is ~ 0.05 around 2 = 1. Therefore, if we 
want to study the evolution of the cosmic acceleration in 
a precious manner, we need to have a resolution with an 
accuracy level AR ~ 0.005, assuming a signal-to-noise 
ratio of 10, around a redshift interval dz ~ 0.1. With 
this sensitivity for the redshift change, Az, we can set 
the upper limit for Q.gw as 

n GW ~ 4.0 x 10 2 x 10" 1 x (5 x 10~ 3 ) 2 - 10~ 3 , (17) 

where the factor 10 _1 represents an average estimate on 
the improvement due to the increase in the effective num- 
ber of the shells. If shells between redshift of 1 and 3 can 
be divided at intervals of 0.1, this increase is a factor of 
a few larger and one further improves the limit on Hqw- 

III. COSMOLOGICAL PROBES 
A. Redshift variations in the line emission 

The use of time variation associated with the redshift 
of distant sources to measure global cosmological param- 
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eters is described in Ref. 19]. The resulting change is an 
overall common shift to the redshift distribution and the 
resulting change can be obtained by data separated in 
time. For this purpose, the most useful observations are 
the ones that involve narrow lines either in emission or 
absorption since the precise redshift can easily be ascer- 
tain from spectroscopic observations. In Ref. [l9| . the use 
of the Lyman-a forest in absorption towards luminous 
quasars has been discussed. The expected mean shift in 
the Lyman-a forest, at a redshift ~ 3, is roughly 2 m s -1 
over a time period of 100 years, and for presently favor- 
able cosmological parameters. While this is a small num- 
ber, spectroscopic variations at the level of ~ 1 m s _1 are 
now routinely measured with 10 meter class telescopes to- 
wards bright stars in extra-Solar planetary searches using 
radial velocity measurements [27| . The advent of 30 me- 
ter class telescopes over the next decade or more and 
the resulting increase in spectral resolution will improve 
searches for small velocity shifts in high-redshift quasar 
spectra. 

While the precise measurement of the overall cosmo- 
logical shift in a single quasar will require observations 
that span many tens of years, a detection of the global 
change can be obtained through a statistical study of 
a large number of quasars over a short time interval. 
While this will allow a measurement of the cosmological 
parameters |T^ . the effect due to a gravitational wave 
background is high order and requires a study related to 
the quadrupolar pattern of redshift shift. While the use 
of a quadrupolar pattern allows one to separate the ef- 
fect related to cosmological acceleration, which is present 
in the monopole, and that due to local motion, present 
in the dipole, with that expected from a background of 
gravitational waves, the measurement, especially for an 
expected low energy density for the gravitational wave 
background is harder. For Hqw at the level of 10~ 10 
the anisotropy is at the level of 10 -5 to 10~ 6 , so the di- 
rect constrain one can put on Oqw depends how well the 
monopole can be established. 

Using estimates in Ref. for observations of 1000 
quasar sight lines with improved spectroscopic data with 
a pixel scale of 0.5 km sec -1 ) over a 10 year span, and 
assuming no systematic or additional statistical errors 
such as in the calibration of the spectrum, we find that 
one can constrain the monopole to an accuracy of ~ 1%. 
The resulting constrain on the low frequency, with wave- 
lengths out to z ~ 3, background of gravitational waves 
is at the level ilcw < 10~ 3 . While the constrain is not 
significant, it is still useful to obtain this directly from 
the data as this will be one of the few methods to con- 
strain the energy density of gravitational waves around 
10" 12 Hz. 

In addition to the Lyman-a forest, one can constrain 
the energy density of gravitational waves using the red- 
shift distribution of 21 cm emitter. Here, the line emis- 
sion is at low radio wavelengths and statistics are ex- 
pected to significantly improve given the expected narrow 
width of these lines, especially before the era of reioniza- 



tion, and the narrow band observations one can, in princi- 
ple, achieve at low radio frequencies. Expected statistics 
of HI emitters before reionization is not well known and 
the complication here will be more related to foreground 
contamination 28] more than the number of objects one 
can use to make these measurements. Improvements in 
foreground removal techniques suggest that the contam- 
ination can be mostly remove d l29j . 

We follow estimates in Ref. _30j which assume a model 
for the cosmological distribution of neutral gas before 
reionization based on the halo approach ,31] and using 
so-called "mini halos" 32] that are beginning to establish 
with adiabatically cooled gas with a temperature below 
that of the CMB thermal radiation at those redshifts. We 
estimate the line width of each signal to be around 3 km 
sec -1 and to be dominated by thermal motions within 
each halo. This compares to line widths of order 20 km 
sec -1 in the Lyman-a spectrum. Since the characteristic 
mass scale of each of these mini halos is around 10 5 M sun 
to 10 6 M sun , one expects a total of, in principle, 10 18 
halos across the whole sky with a size of 1 kpc, or a 
projected size of 30 milliarcseconds. While the whole 
set of mini halos are not needed, using a sample of 10 6 
brightest 21 cm lines over the whole sky, and assuming 
observations with an observational bandwidth of 0.2 kHz, 
which corresponds to a velocity of ~ 0.4 km sec -1 , we find 
that one can, in principle, constrain the background to 
be below SIgw < 10~ 5 . The ultimate limit related to 21 
cm anisotropics, using the whole sample of 10 18 halos and 
assuming sufficient resolution to resolve them in redshift 
space and ignoring source confusion, is Qgw < 10 
with lOyr observation. This is comparable to the current 
limit from CMB anisotropics, but at frequencies between 
10~ 18 Hz and 10~ 16 Hz. 



B. Gravitational waves from neutron star binaries 
around 1Hz 

Here we discuss the prospect of the measurement of 
the background Qgw at the very low frequency regime 
(< 10 _12 Hz) with proposed space gravitational wave mis- 
sions, such as, the Big Bang Observer (BBO) or DE- 
CIGO, whose optimum sensitivities are between LISA 
and ground based detectors and around 0.1 Hz to 1 Hz. 
Note that this frequency range of 0.1 Hz to 1 Hz is largely 
different from the band at / < 10~ 12 Hz. 

The most important aim of these detectors is the di- 
rect detection of the gravitational wave background from 
early universe with Qgw ^ 10~ 15 level at 0.1 to 1 Hz. For 
this measurement, the foreground gravitational waves 
from various astrophysical sources must be subtracted 
from the data streams down to some appropriate level. 
The cosmological neutron star binaries are expected to 
produce a strong foreground, and we need an extensive 
fitting analysis with accurate templates to remove these 
sources from the data. In this fitting procedure we can 
measure not only the intrinsic binary parameters, but 
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also the total extrinsic effect made by the cosmic or the 
local acceleration through the phase variation related to 
X. In the presence of a low frequency gravitational wave 
background, extracted X values for a sample of bina- 
ries spread over the sky, over a certain bin width in 
redshift, are suppose to vary with an anisotropic pat- 
tern. As discussed, the quadrupolar pattern can in re- 
turn be used as a probe of the low frequency gravita- 
tional wave background. Therefore cosmological neutron 
star binaries can be regarded as probes of the dark en- 
ergy, th rough the monopole variation related to X(z) 
[22. l23| , or the low frequency gravitational wave back- 
ground through the anisotropy pattern of X. While the 
discussion is related to neutron star binaries detectable 
in the 0.1 Hz to 1 Hz band, we can also make a simi- 
lar argument for stellar mass black hole binaries in the 
same band, though their intrinsic waveforms would be 
generally more complicated and the extraction of X(z), 
through phase changes, would be a complicated proce- 
dure. The local acceleration might be generally larger 
than the neutron star binaries. 

At present, it is not clear how well we can actually re- 
solve and subtract the gravitational wave contribution of 
each merging neutron star binary in BBO or DECIGO 
data streams. For example, the confusion effect is more 
important for a binary with a longer time before coales- 
cence, as its gravitational waves are in a lower frequency 
regime where the number density of binaries per unit 
frequency increases. In addition to binaries, the back- 
ground made by other sources such as supernovae could 
be a fundamental noise |33| . 

Now we study how well we might set the constraint for 
flaw at / < 10~ 12 Hz based on the recent analysis for the 
measurement of the cosmic acceleration R = X(z)/Hq 
|23j. Following their paper, we use the coalescence rate 
of neutron star binaries at 10 _6 /yr/Mpc 3 assuming it 
is constant in time. The total number of binaries in 
a redshift shell z ~ z + dz is ~ 10 4 /yr for a width 
dz ~ 0.1 at z ~ 1. As we mentioned, in addition to 
the measurement error AXdet caused by detector noises, 
the random local acceleration becomes an effective noise 
AXi oca i for estimating the cosmic acceleration signal 
or the gravitational wave background. In the case of 
a neutron star binary, its magnitude would correspond 
to AXinrai/Ho ~ 10 _1 in terms of the noise AR for 
each binary |3J|. Therefore, for a shell with a width 
dz ~ 0.1, the estimation of the common signal is limited 
by the fluctuations d ue to the local acceleration around 
(Aii)sheU ~ 10 _1 x Vl0~ 4 ~ 10 -3 for an observation pe- 
riod of order of one year. The final fluctuations after us- 
ing all shells becomes AR~ (Ai?) s hellXV7Vshell ~ 10~ 3 ' 5 . 
With an effective number for shells of A^eii = 10, we 
can set a upper limit of flow ~5x 10~ 5 using the rela- 
tion ljPo|) . 

As for the measurement error AXdet determined by 
the detector noises, we can reach the level (Ai?) s h c ii ~ 
3 x 10 -4 (corresponding to AXdet/ Hq ~ 3 x 10 -2 for each 
binary) with a 3 year integration by using the ultimate 



DECIGO (see figures 3 and 4 in Ref. (H). But with 
a one year integration, the measurement error becomes 
(Ai?) s hoii ~3x 10 -3 and dominates the fluctuation by 
the local acceleration. In this estimation we did not in- 
clude astrophysical confusion noises. 



C. Small Angular Scale CMB Anisotropies 

If a gravitational wave background exist at frequencies, 
say, between 10 -14 Hz and 10~ 12 Hz, the time evolution 
of the background would produce a signature in CMB 
anisotropies through metric perturbations |35j . This is 
similar to temperature fluctuations produced at large an- 
gular scales with horizon-scale waves. In terms of the 
angular power spectrum of CMB anisotropies, the signa- 
ture would be at multipole moments between ~ 6 x 10 4 
and 6 x 10 6 . Such small-scale anisotropies, at arcsec- 
ond angular scales, are not probed by any of the present 
observations. In fact, no dedicated plans exist for the 
CMB anisotropy observations at such small scales though 
certain radio interferometers planned for purposed other 
than CMB, such as the Atacama Large Millimeter Tele- 
scope (ALMA), can be used to make the required mea- 
surements. 

The fluctuations at these scales are dominated by un- 
resolved radio point sources, which produce a shot-noise 
spectrum. For comparison, assuming a radio point- 
source removal similar of the Cosmic Background Imager 
(CBI; (H3 } experiment, where residual fluctuations were 
at the level of 150 /iK 2 in the smallest angular scale bin in 
multipole moments between 2000 and 3500, we find that 
the resulting limit to CMB temperature fluctuations in a 
bin with multipole moments between ~ 6x 10 4 and 6x 10 6 
is AT < 3000/iK. This results in a limit for the gravita- 
tional wave background between 10~ 14 Hz and 10 -12 Hz 
of Ogw < 10~ 9 . To reach the ultimate level implied by 
21 cm observations require a factor of 100 improvement 
in foreground removal. Since foregrounds are already re- 
stricted to a few percent level, this is a challenging task. 



IV. SUMMARY 

The curl modes of cosmic microwave background po- 
larization allow one to indirectly constrain the primor- 
dial background of gravitational waves with wavelengths 
roughly the horizon size or larger with frequencies below 
10~ 16 Hz. The planned high precision timing observa- 
tions of a large sample of millisecond pulsars with the 
Pulsar Timing Array or with the Square Kilometer Array 
can either detect or constrain the stochastic gravitational 
wave background at frequencies greater than roughly 0.1 
years -1 . While there are no strong observational con- 
straints on the gravitational wave background across six 
or more orders of magnitude between 10~ 16 Hz and 10~ 10 
Hz, we suggest that by monitoring the anisotropy pat- 
tern of time variation of the redshift related to a sample 
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of high redshift objects one can constrain the gravita- 
tional wave background below 10~ 12 Hz where methods 
using Galactic objects do not work well. Direct mea- 
surement of the time variation of redshift is one of the 
potential method to get information of the dark energy. 
Analyses in this paper would be helpful to discuss the 
prospects of the low frequency gravitational wave back- 
ground measurement in a given observational condition 
of the dark energy study. Useful observations for the 
monitoring of an anisotropy signal in a global redshift 
change include spectroscopic observations of the Ly-a 
forest in absorption towards a sample of quasars, red- 
shifted 21 cm line observations, either in absorption or 



emission, towards a sample of neutral HI regions before 
or during reionization, and high frequency (0.1 Hz to 1 
Hz) gravitational wave analysis of a sample of neutron 
star — neutron star binaries detected with gravitational 
wave instruments such as the Decihertz Interferometer 
Gravitational Wave Observatory (DECIGO). For reason- 
able observations in the future, we find best limits at the 
level of Hgw < 10~ 5 at a frequency around 10~ 12 Hz. 
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